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Abstract 


A  system  was  developed  which  could  measure  the  instantaneous 
velocities  in  the  wake  of  a  model  helicopter  rotor  operatinp  at  low 
advance  ratios.  A  three-wire  hot  film  probe  was  mounted  on  a  traverse 
and  placed  at  many  positions  in  the  wake.  The  output  of  the  probe  was 
fed  into  an  on-line  computer  operating  in  an  interactive  mode.  Computer 
generated  vector  plots  were  made  of  both  time  averaged  velocities  and 
instantaneous  velocities  for  the  case  of  y  =  0.06. 


CHAPTER  I 


INTRODUCTION 


Recent  requirements  for  an  Advanced  Attack  Helicopter 
(AAH)  whose  primary  missions  will  be  anti-armor  in  a 
battle  zone,  and  low  speed  terrain  flight  (Nap-Of-The- 
Earth,  NOE)  has  increased  the  need  for  data  of  helicopter 
rotor  wakes  at  low  advance  ratios. 

The  blade  tip  vortices  constitute  the  primary  distur¬ 
bances  in  the  rotor  wake,  and  contribute  significantly  to 
the  dynamic  loads  and  acoustic  characteristics  of  the 
helicopter.  Varitvs  sophisticated  theories  have  been 
developed  to  calculate  the  velocities  of  the  wake,  Refs.  (1-6) 
the  vortex  positions,  and  the  resulting  dynamic  loads. 

These  have  been  compared  to  measured  loads,  but  it  is  not 
known  whether  the  discrepancies  are  due  to  nonlinear  and 
3-dimensional  effects  or  to  inaccuracies  in  the  predicted 
vortex  characteristics.  It  is,  therefore,  desirable  to 
obtain  accurate  measurements  of  the  rotor  wake,  the  vortex 
position,  and  vortex  characteristics. 

Early  measurements  of  the  flow  field  used  total 
pressure  tubes  (Ref.  7).  These  results  were  very  useful 
in  demonstrating  the  nonuniformity  of  the  time  averaged 
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flow.  However  the  instantaneous  velocities,  which  are 
especially  important  for  accurate  loads  calculations,  were 
not  measured.  Hot-film  anemometers  have  been  used  for  rotor 
wake  studies  and  reported  in  Refs.  8,  9,  and  10.  The 
results  of  Ref.  8  allowed  the  structure  of  a  tip  vortex  to 
be  defined  for  the  hover  conditions  tested,  whereas  in 
Ref.  9  and  10  the  averaged  velocities  at  low  advance  ratios 
were  measured.  Photographic  techniques  have  also  been  used. 
Ref.  11  and  12  show  results  from  two  investigations  using 
Schlieren  photography  systems.  Water  and  smoke  studies 
using  photography  are  presented  in  Ref.  13  and  14.  Both  of 
these  methods  are  useful  for  flow  visualization,  but  data 
reduction  is  somewhat  difficult  and  accuracy  is  not  high. 

The  first  demonstrated  use  of  a  Laser  Doppler  Anemometer 
(LDA)  for  rotor  investigation  was  that  in  Ref.  15.  In  more 
recent  works  with  an  LDA,  such  as  Ref.  16,  instantaneous  as 
well  as  time  averaged  velocities  were  measured.  The  laser 
system  allows  flow  measurements  to  be  made  very  close  to 
the  blade,  and  the  vortex  velocity  distribution  and  core 
size  were  determined. 

Very  little  is  known  about  the  transition  flight  speed 
(0  <  y  <  0.2)  performance.  At  these  speeds,  very  large 
wake  distortion  and  very  severe  vibratory  loads  occur.  The 
purpose  of  this  work  was  to  develop  the  instrumentation  and 
the  experimental  methods  needed  to  measure  the  instantane¬ 
ous  velocities  in  the  vicinity  of  a  small  helicopter  rotor 
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in  a  wind  tunnel  operated  at  low  advance  ratios  and  to  make 
measurements  of  typical  flow  fields.  The  measurements  were 
made  with  a  three  dimensional  hot-film  anemometer.  The 
data  acquisition  was  done  by  a  DEC  PDP  11/60  minicomputer 
with  a  GENRAD  Analog  to  Digital  System.  A  three  dimensional 
traverse  system  was  used  to  traverse  the  anemometer  probe 
throughout  the  rotor  wake. 


CHAPTER  II 


BACKGROUND 

One  can  distinguish  between  two  different  types  of 
velocity  information  at  any  point  in  the  rotor  wake.  The 
first  type  is  the  time-average  velocities  induced  by  the 
helicopter's  rotor  in  horizontal  flight,  which  is  the  con¬ 
stant  term  of  the  Fourier  representation  of  the  wake.  The 
second  type  of  velocity  information  is  the  instantaneous 
velocities  induced  by  the  rotor. 

In  a  steady-state  region  of  flight,  the  instantaneous 
velocities  are  periodic  with  respect  to  time.  The  period 
is  equal  to  the  time  required  for  the  rotor  to  turn  through 
an  angle  equal  to  that  between  two  consecutive  blades.  In 
the  consideration  of  flow  around  aerodynamic  surfaces, 
fixed  with  respect  to  the  helicopter,  such  as  a  wing, 
stabilizer,  or  fuselage,  instantaneous  induced  velocities 
are  usually  replaced  by  the  time-average  values  at  a  given 
point.  The  smaller  the  pulsations  are  in  comparison  with 
their  time-average  values,  the  more  correct  this  hypothesis 
would  be. 

The  theoretical  method  of  predicting  the  instantaneous 
induced  velocities  is  based  on  the  fundamental  assumption 
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of  the  linear  theory,  according  to  which  the  transport  of 
the  vortex  elements  occurs  with  constant  velocity  (Ref.  17) . 
The  linear  theory  gives  good  results  as  long  as  the  induced 
velocities  are  small  in  comparison  with  those  of  the  undis¬ 
turbed  flow  approaching  the  rotor.  In  that  case  the  velo¬ 
city  of  the  transportation  of  vortices  is  not  much  different 
from  that  of  the  undisturbed  flow.  This  case  is  typical  for 
high  speed  helicopter  flight. 

No  theoretical  treatment  is  presently  available  for 
low  advance  ratios,  because  the  phenomenon  studied  actually 
deviates  significantly  from  the  linear  model.  For  this 
case  the  magnitudes  of  the  induced  velocities  are  comparable 
to  those  of  the  flight  itself  and  large  downwash  angles 
occur.  Further  complicating  the  problem  in  the  low-speed 
flight  region  are  the  numerous  blade-vortex  interactions 
that  occur  (Ref.  18) .  In  order  to  improve  the  prediction 
accuracy  of  the  rotor's  wake  in  those  regimes  of  flight,  it 
is  desirable  to  have  a  better  insight  in  the  actual  physi¬ 
cal  phenomena  that  occurs.  This  work  is  an  attempt  to 
provide  the  means  of  producing  data  for  helicopter  rotors 
in  forward  flight  at  low  advance  ratios. 


CHAPTER  III 


EXPERIMENTAL  APPARATUS  AND  PROCEDURE 


o.l  The  Test  Rig 

The  tests  were  conducted  in  a  wind  tunnel  with  a  test 
section  of  8  ft  x  4  ft,  in  The  Department  of  Mechanical 
Engineering  at  The  Ohio  State  University.  The  flow  turbu¬ 
lence  is  reduced  by  means  of  2  to  1  contraction  ratio  and 
nested  1/8  inch  diameter  straws.  The  wind  tunnel  air 


velocity  is  variable  up  to 
The  tests  were  carried 
rotor  with  twisted  blades. 
Table  3.1 

Table  3.1  Detai 
Number  of  Blades 
Rotor  Radius,  ft 
Blade  Chord,  in 
Rotor  Solidity  bc/rrR 
Root  Cutout  %  R 
Blade  Taper  Ratio 
Coning  Angle 
Blade  Twist 


5  ft/sec. 

out  on  a  teetering  two  bladed 
Details  of  the  rotor  are  in 

Is  of  the  Rotor 
2 

1.25 

2.1875 

0.0928 

12.08 

1 

0° 

8°  (From  12.08%  R  to  100%  R) 
0.004260  slug-ft2 
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gure  3-2  Dimensions  of  the  Test  Section 


Airfoil  Section 


NACA  0012 


The  rotor  is  driven  by  a  variable  speed  electric  motor 
The  speed  is  variable  up  to  a  maximum  of  3000  RPM  which 
results  in  tip  speeds  up  to  400  ft/sec.  The  collective 
pitch  and  the  rotor  shaft  tilt  angle  are  variable  but  are 
constant  for  each  case  run. 

3.2  Data  Acquisition  System 

3.2.1  General  Description 

The  data  acquisition  system  consists  of  two 
basic  subsystems:  1.  the  velocity  measurement  system  and 
2.  the  minicomputer  and  the  A/D  system. 

The  velocity  measurement  system  is  located  at  the 
windtunnel  on  the  first  floor  of  The  Mechanical  Engineering 
building  and  consists  of  a  three-dimensional  hot-film  probe 
connected  to  three  constant  temperature  anemometer  circuits 
The  probe  is  attached  to  a  traverse  mechanism  positioned 
in  the  wind  tunnel  enabling  the  operator  to  locate  the 
probe  at  the  desired  point  in  the  tested  volume. 

The  minicomputer  and  the  A/D  are  located  in  the 
computer  room  (ADML  laboratory)  on  the  second  floor  of  The 
Mechanical  Engineering  building.  The  operator  interfaces 
the  computer  by  means  of  a  alphanumeric  terminal.  The 
triggering  of  the  A/D  is  done  from  the  wind  tunnel  by  a 
special  interface  designed  and  built  for  this  purpose. 


Computer  Room 
2nd  Floor 
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Figure  3-4  Block  Diagram  of  the  Anemometry 


Figure  3.3  shows  the  interconnection  of  the  entire  system. 

3.2.2  The  Velocity  Measurement  System 

The  rotor  induced  velocities  were  measured  by 
a  three  dimensional  hot-film  probe  (TSI  1294-20-18)  connect¬ 
ed  to  three  DISA  55D05/102C  constant  temperature  anemometer 
circuits  w^th  a  SOLA  83-15-2216  power  supply.  In  order  to 
improve  the  frequency  response  of  the  hot -film  anemometer 
a  lead  compensator  was  added  to  each  channel  (Appendix  A). 
Because  of  the  mode  of  operation  of  the  A/D  converter 
channels  B  and  C  had  to  be  delayed,  6  y  sec  and  12  y  sec 
respectively,  in  order  that  the  computer  will  "see"  the 
same  angle  of  the  rotor  in  time.  The  delay  lines  used  were 
Allen  Avionics  1000  ohms  impedance  6  y  sec  B06P0Z01K.  One 
delay  line  was  used  for  ' B ’  and  two  6  y  sec  delay  lines 
in  series  were  used  for  ’C  to  create  the  12  y  sec  delay. 

The  output  from  the  anemometer  thru  the  delay  lines  was 
connected  to  three  DANA  2210  DC  amplifiers.  Figure  3.4 
shows  the  complete  velocity  measuring  system  including 
the  lead  compensators  and  the  delay  lines. 

3. 2. 2.1  Directional  Sensitivity  of  Cylindrical 
Sensors 

The  basic  hot-wire  hot-film  anemometer 
output  is  a  voltage  which  is  related  to  the  fluid  velocity 
approximately  as: 


I 
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E2  «  A  ♦  B  U1/2  (3-1) 

where  U  is  the  stream  velocity,  E  is  the  voltage  output 
and  A,  B  are  constants  determined  by  the  electronic  circuit, 
the  physical  properties  of  the  sensor  and  by  the  fluid 
properties . 


If  the  fluid  velocity  is  not  perpendicular  to  the 
sensor,  equation  (3-1)  does  not  hold,  since  the  heat 
transfer  from  the  sensor  is  reduced.  A  convenient  way  to 
overcome  this  problem  is  to  define  an  "effective  cooling 
velocity"  as  the  perpendicular  velocity  which  produces  the 
same  cooling  effect  as  the  actual,  non -perpendicular 
velocity  and  to  correlate  it  to  the  true  stream  velocity. 
The  simplest  relationship  suggested  was: 


U  cos  0 


(3-2) 


13 


Equation  (3-2)  neglects  the  contribution  of  the  flow  along 
the  sensor  and  is  not  accurate  for  $  >  40w.  More  accurate 
relations  have  been  suggested  in  References  22,  23,  24. 

For  this  work  the  relationship  proposed  by  Ref.  (23)  was 
used  after  an  extensive  investigation  in  Ref.  21. 

This  relationship  is  defined  by  eq.  3-3 

Veff  =  u2(cos2  *  +  k2  sin2*)  (3-3) 

where : 

V  ff  is  the  effective  cooling  velocity  normal  to  the 
sensor . 

U  is  the  stream  velocity. 

0  is  the  angle  between  the  free  stream  and  the  normal 
to  the  sensor. 

k  is  a  dimensionless  constant. 

3. 2. 2. 2  Data  Reduction  Method  for  the 
3-Dimensional  Probe 

The  three  dimensional  probe  consists 
of  three  orthogonal,  cylindrical  sensors.  (Equation  3-3 
holds  for  each  sensor.)  The  effective  cooling  velocities 
of  the  three  sensors  are: 

V2  =  U2  (cos 2$i  +  k2sin24>i)  i  -  1,2,3  (3-4) 

Adding  the  three  above  equations  leads  to: 

ZV2  *  U2 (2  ♦  k2) 


(3-5) 


I 


Sensor  C 


Figure 


(b) 

3-6  Geometry  of  a  3-D  Hot-Film  Probe 


(a)  Sensors 


(b)  Velocities 
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where  the  following  trigonometric  relationships  were  used: 


Icos 


Isin2<J)i  *  1 


(3-6) 

(3-7) 


Solving  for  the  mean  velocity: 


.  ivf 
u  - ^T 


(3-8) 


substituting  equations  (3-8)  in  equations  (3-4)  the  direc¬ 
tional  angles  can  be  obtained 


s  in.  0 


Vi 

i  (I  *  I~1  1/(1  *  k2)  i  =  1,2,3 


(3-9) 


equations  (3-8)  and  (3-9)  precisely  define  the  total  velo¬ 
city  vector  and  its  direction  in  the  probe  system  of 
coordinates . 

2 

The  above  analysis  is  correct,  provided  that  k  is  the 
same  for  all  three  wires  and  that  it  is  not  a  function  of 
the  velocity  vector  U.  fig-  3.7  shows  a  typical 
relationship  between  k  and  U.  It  is  apparent  that  at 
higher  velocities  k  can  be  assumed  constant  whereas 
at  lower  velocities  this  assumption  will  lead  to 
erroneous  results. 

To  overcome  the  above  stated  problem  an  iterative 

2 

program  was  developed  to  account  for  the  dependency  of  k 
on  the  total  velocity  vector  U.  For  this  case 
the  following  analysis  was  used 
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Extensive  investigation  on  the  influence  of  k  on  U  and  <f>^ 

2 

showed  that  k  had  a  negligible  effect  on  the  calculations 

of  U,  and  a  significant  influence  on  the  angles.  Thus,  to 

simplify  the  iterations  a  constant  value  of  k  (.based  on 

2  2  7 

the  average  values  of  k^,  kg,  k£  was  used  in  calculating 

2 

the  total  vector,  and  a  local  value  of  was  used  in 

2 

calculating  <J>^  .  The  determination  of  k  =  f(U)  is  explained 
in  Section  3. 2. 2. 4. 2  and  the  iterative  computer  program  is 
shown  in  Appendix  B. 

3. 2. 2. 3  Transformation  of  Velocities  from  the 
Probe  System  of  Coordinates  to  the 
Tunnel  System  of  Coordinates 
The  probe  used  was  a  T.S.I.  1294-2U-18 
hot-film  probe  shown  in  Figure  3.8.  The  three  sensors  are 
mutually  perpendicular  and  inclined  54.74°  to  the  probe 
axis.  The  spatial  orientation  of  the  probe  in  the  wind 
tunnel  is  shown  in  Fig.  3.9,  and  the  probe  coordinate 
system  is  shown  in  Fig.  3 . 10 ,  relative  to  the  rotor  coordi¬ 
nate  system.  The  orientation  of  the  probe  is  such  that 
axis  C  lies  in  plane  X  4  Z  and  plane  BOAD  is  perpen¬ 
dicular  to  the  same  plane.  For  this  orientation  we  can 
perform  the  coordinate  transformation  knowing  that: 

<EOC  =  35.26“ 

<E0D  =  54.74° 


<AOD  =  <BOD  =  45.0° 


Detail  \  support. 


3  Mutually 

Perpendicular  Sensors 
Symmetrically  Positioned 
with  Respect  to  Probe 
Axis 


Figure  3-8  TSI  1294-20-18  Three-Dimensional  Hot-Film 
Probe 
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Figure  3-9  The  Spatial  Orientation  of  the  Probe  in  the 
Wind  Tunnel 


Figure  3-10  Probe  Coordinate  System  Relative  to  Wind 
Tunnel  Coordinate  System 
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Equations  3-13  through  3-16  define  this  transformation. 

4  V  -  4-  0.57735(11.  /4  (3-13)  '  X. 

7  V/=  0.7071  (U  -  U  ).  v/  -  4  (3-14) 

/ 

V  =-0.4082(U.  +  Ufa)  -  0.816bUc  V  (3-15) 

'x  * 

3. 2. 2. 4  Calibration  of  the  Three-Dimensional 
Probe 


The  calibration  of  the  three-dimension - 
al  probe  consists  of  two  stages.  The  first  is  to  obtain  /  -•( 

?  Corf1*' 

the  voltage  versus  velocity  curve  for  each  sensor  and  the  l 
second  is  to  determine  the  dependency  of  the  yaw  angle  J 


sensitivity  constant  k  ,  on  the  total  vector  U. 


3. 2, 2. 4.1  Thf*  Velocity  Calibration 

The  velocity  calibration  is 

performed  on  a  calibrator  which  consists  of  a  nozzle  with 
exit  diameter  of  U.75  inches,  a  turbulence  reduction  cham¬ 
ber  and  a  probe  holder  that  can  vary  the  angle  of  the  probe 
relative  to  the  air  stream.  The  air  is  supplied  from  a  100 
P.S.I.  supply  line  and  controlled  by  a  needle  valve.  An 
investigation  done  on  this  calibrator  in  Ref.  21  showed 
that  0.4  inches  from  the  exit  plane  the  jet  had  a  uniform 
velocity  core  of  0.6  inches.  This  location  was  used  to  get 
a  correlation  curve  between  the  nozzle's  velocity  and  the 
chamber  static  pressure.  (Pig.  3.11).  The  same  location 
was  used  for  the  hot-film  anemometer  calibration. 


Jet  Velocity  @  78°F  -  ft/sec 


0.0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1.0 


The  Chamber  Static  Pressure  -  inches  1^0 

Figure  3~11  The  Calibration  Curve  for  the  Velocity 
Calibrator 

The  hot  film  anemometer  was  calibrated  by  adjusting 
the  chamber's  static  pressure,  which  corresponds  to  a 
specific  velocity,  and  recording  the  anemometer's  output 
voltage.  This  calibration  was  done  in  the  range  of  0  to  70 
ft/sec,. 

The  author  in  Ref.  21  has  shown  that  the  sensors' 
prongs  do  effect  the  "effective  cooling  velocity".  This 
velocity  is  1  to  3  ft/sec  smaller  when  the  prongs  are  in 
the  plane  of  the  velocity  vector  than  when  the  prongs  are 
perpendicular  to  that  plane,  for  the  range  of  10  to  100 
ft/sec.  For  the  probe  orientation  in  the  wind  tunnel, 
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sensors  A  and  B  had  their  prongs  perpendicular  to  the  total 
vector,  most  of  the  time,  and  sensor  C  had  its  prongs  in 
the  plane  of  the  velocity  vector.  Therefore,  sensors  A 
and  B  were  calibrated  with  their  prongs  perpendicular  to 
the  jet  velocity  and  sensor  C  was  calibrated  with  its 
prongs  in  the  plane  of  the  jet.  The  resulting  calibration 
curves  are  shown  in  Fig.  3.12. 

Since  the  calculations  of  the  velocity  was  done  by  a 
computer,  a  forth  order  polynomial  was  fitted  to  the  cali¬ 
bration  curves  of  each  sensor.  The  effective  velocities  were 
calculated  with  this  equation,  as  shown  in  eq.  3-16. 

Veff  *  C4e4  *  C3e3  *  C2e2  *  C1E  *  C0  '3-16) 

The  constants  of  the  polynomials  were  calculated  by  an 
existing  program  using  the  Method  of  Least  Squares.  The 
curve  fit  matched  the  experimental  results  within  0.4 
ft/sec . 

3. 2. 2. 4. 2  Directional  Sensitivity 
Calibration 

The  purpose  of  this  cali- 
2 

bration  was  to  determine  the  dependency  of  k  on  the 
velocity  vector.  An  extensive  investigation  done  for  this 
purpose  revealed  that  k  is  strongly  dependent  on  the 
wire's  Reynolds  number  where 


Re 


yd 

v 


(3-17) 


siioa  ‘  3DvnoA  indino 


Calibration  Curves 
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However,  d  and  v  are  constant  for  each  case,  and  Rg  becomes 

7 

a  function  of  U  alone.  Therefore,  k"  was  evaluated  also  as 
a  function  of  U.  For  high  Reynolds  numbers,  k  is  nearly  a 
constant.  Thus,  only  by  changing  wire  diameter  (different 
probes)  the  values  of  k2  would  change  even  though  the  regime 
of  velocity  was  not  changed.  For  the  TSI  1294-60-18  probe 
with  U.006  inches  wire  diameter,  a  single  value  of  k  was 
adequate,  similar  to  the  method  discussed  in  Ref.  21. 

However,  for  the  measurements  of  instantaneous  veloci¬ 
ties,  a  finer  wire  was  used  in  order  to  improve  the  system 

frequency  response.  For  this  purpose  a  TSI  1294-20-18  probe 

2 

was  used  with  0.002  inches  wire  diameter.  For  this  probe  k 

showed  a  strong  dependency  on  the  velocity  (lower  Reynolds 

number  due  to  smaller  diameter).  As  seen  in  Figure  5-7. 

2 

In  order  to  determine  the  correlation  between  k  and  U, 
the  following  procedure  was  used.  Rewritemg  equation  (3-3) 
as : 

V  2 

(-^— )  «  cos2<j>  +  k2sin2<j>  (3-17) 

2 

one  can  obtain  the  values  of  k4-  by  calculating  the  coeffi¬ 
cients  of  equation  (3-17)  when  the  ratio  is  given 

for  a  specific  <J> ,  This  was  done  by  holding  U  constant  and 
changing  4  from  0  to  90  degrees.  Repeating  this  procedure 
for  U  in  the  range  of  10  to  70  ft/sec.  the  dependency  of 

k2  on  U  was  obtained.  Then  the  data  obtained  was  curve 

2 

fitted  so  that  the  relationship  k  =  f(U)  could  be  used  m 
the  computer  program.  Typical  correlations  are  shown  in 
equations  3-18  through  3-20. 


For  sensor  A: 
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k2A  =  2.546x10" 7xU4  -4 . 533xl0"5xU3+0 . 00  3xU2 - . 08955xU+l . 104  2 

(3-18) 

For  sensor  B: 

k2„ =  1. 0 33x10 "7xU4-l . 8787xl0~3xU3+. 001296 xU2-.U41xU+. 54244 

(3-19) 

For  sensor  C: 

k2c  =  2. 4781xl0'7xU4-4. Z8  28xl0'5xU3+. 002  75  5xU2- .07928XU1- . 9092 

(3-20) 

The  incorporation  of  tnis  relationship  with  the  complete 
computer  program  and  the  iterative  computation  are  shown 
in  Appendix  B. 

The  authors  of  References  23  and  24  have  found  that 

2 

k  is  dependent  of  the  geometry  of  the  sensor  (aspect 
ratio  1/d)  and  almost  independent  of  other  properties  of 
the  sensor  and  the  fluid.  Therefore,  the  value  of  k2  is 
not  expected  to  change  with  time  due  to  contamination  and 

corrosion.  Intuitively,  a  similar  assumption  can  be  made 

2 

about  the  dependency  of  k  on  U,  thus  the  calibration 
curves  obtained  in  this  work  for  k2  *  f(U)  may  be  used  in 
future  applications  without  recaiibration  and  it  is  good 
for  the  range  of  velocities  of  10  to  70  ft/sec. 

3.2.3  The  Minicomputer  and  the  A/D  System 

The  new  minicomputer,  installed  in  The  Depart¬ 
ment  of  Mechanical  Engineering  at  The  Ohio  State  University, 
in  the  summer  of  1978,  provided  the  means  of  improving  the 
system  of  average  data  acquisition  and  extending  the  work 
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to  instantaneous  velocity  measurement. 

3.2. 3.1  The  Minicomputer 

The  minicomputer  is  a  Digital  Equipment 
Corporation  (DEC)  PDP  11/60  and  is  well  suited  to  Fortran 
based  calculations.  It  has  a  2b6  KBYTES  memory  (RAM)  and 
64  KBYTES  are  available  to  the  user.  The  PDP  hosts  gra¬ 
phics  (Tektronix  4014-1)  and  alphanumeric  (Infoton)  termi¬ 
nals.  The  terminals  are  supported  by  a  high  speed  Print- 
ronix  Printer -Plotter  and  a  Tektronix  hardcopy  unit.  Data 
and  program  storage  are  provided  by  a  dual  floppy  disc 
drive  and  two  cartridge  disc  drives. 

3. 2.3. 2  The  Analog  to  Digital  System  (ADS) 

Two  A/D  converters  became  available 

during  the  program.  The  first  was  the  DEC  LPA  system  which 
can  sample  up  to  40  KHZ  and  can  be  multiplexed  up  to  16 
channels.  The  second  was  the  GENRAU  4  channel  ADS  with 
sampling  frequency  up  to  160  KHZ.  initially  the  DEC  LPA 
system  was  used  because  it  was  planned  to  development 
eventually  automate  the  data  qcquisition.  The  computer 
would  control  the  RPM  of  the  rotor  and  the  x-y-z  location 
of  the  probe.  However,  due  to  hardware  problems  in  the 
LPA  system,  the  data  acquisition  system  was  developed 
around  the  GENRAD  ADS.  (Henceforth  referred  to  as  "the 
ADS") . 
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The  ADS  is  controlled  by  a  canned  program  using  sub¬ 
routines  available  in  the  PDP’s  library.  The  user  has  to 
specify  the  number  of  channels  used  and  the  maximum  expect¬ 
ed  analog  voltage.  This  is  done  by  calling  the  subroutine 
SET ADI ( I POINT ,NCHAN , I RANGE ,KHZ )  wnere:  IPOINT  specifies  the 
number  of  points  to  sample  per  channel.  This  number  is 
encoded  as  shown  in  Table  3.2.  NCHAN  is  the  number  of 
channels  to  sample.  In  our  case  there  are  three  channels. 
IRANGE  is  the  analog  input  voltage  range  encoded  as  shown  in 
Table  3.3,  and  KHZ  specifies  the  corner  trequency  of  the 
antialising  filters  in  the  ADS. 

Table  3-2  Table  3-3 

IPOINT  Code  IRANGE  Code 


CODE  9 

IPOINT 

CODE  # 

ANALOG  RANGE 
(+/-)  VOLTS 

0 

128 

0 

16 

1 

256 

1 

8 

2 

512 

2 

4 

3 

1024 

3 

2 

4 

2048 

4 

1 

5 

4096 

5 

0.5 

6 

8192 

6 

0.25 

7 

16384 

7 

0.125 
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The  subroutine  SETAD1  initiates  the  A/D  conversion,  and  now 
a  second  subroutine  is  called  to  actually  digitize  the 
analog  data.  This  subroutine  is  called  GETAD(BUFFER)  where 
BUFFER  is  an  integer  array  whose  dimensions  are  specified 
by  the  user.  The  digitizing  process  continues  until  the 
array  BUFFER  is  filled.  If  the  triggering  signals  stop 
before  BUFFER  is  filled,  the  computer  will  wait  until  the 
array  is  filled.  Thus,  the  number  of  triggering  signals 
always  has  to  be  larger  than  the  size  of  the  array.  For 
our  case,  it  was  decided  to  average  the  instantaneous  data 
for  15  revolutions  every  5  degrees,  thus  creating  15  x  72  = 
1080  signals.  The  dimension  of  BUFFER  was  selected  to  be 
1024,  which  indicates  that  only  14  complete  revolutions 
will  be  available  for  data  manipulation. 

The  ADS  is  triggered  by  one  pulse  sent  from  the  inter¬ 
face  (see  Appendix  C  for  details) .  The  system  then  samples 
the  first  channel  and  after  6  y  sec  the  next  one,  until  the 
number  of  channels  sampled  equals  the  number  of  channels 
specified  in  the  subroutine  SETAD1.  The  6  y  sec  inter¬ 
channel  interval  caused  significant  errors  in  our  case, 
where  the  three  digitized  values  had  to  correspond  to  the 
same  instant  in  space  and  time.  A  simultaneous  triangular 
wave  of  2000  Hz  had  about  101  error  when  the  digitized 
values  A  and  C  were  compared.  To  overcome  this  problem 
Allen  Avionics  delay  lines  were  used.  Using  these  delay 
lines  the  analog  signal  B  was  delayed  by  6  y  sec,  and  the 
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analog  signal  C  was  delayed  by  12  y  sec.  The  simultaneous 
triangular  wave  test  conducted  with  the  delay  lines  showed 
agreement  between  channels  A  and  C  within  0.26%,  which  is 
of  the  order  of  magnitude  of  the  noise  present  in  the  long 
cable,  connecting  the  ADS  to  the  wind  tunnel. 

The  analog  signals  coming  out  of  the  constant  tempera¬ 
ture  anemometer  were  connected  by  a  S00  ft.  cable  required 
to  get  from  the  wind  tunnel  room,  where  the  experiments  were 
conducted,  to  the  computer  room.  An  analysis,  Appendix  A, 
conducted  on  the  significant  frequency  content  of  the  wake 
at  a  typical  point,  showed  that  the  significant  frequency 
is  contained  in  the  region  below  4000  Hz.  Fortunately,  this 
frequency  was  far  below  the  frequency  where  the  reactive 
elements  of  the  cable  became  significant.  This  frequency  of 
the  cable  was  found  to  be  around  20  KHz.  However,  the  trig¬ 
gering  signal  was  a  square  pulse  of  1  y  sec.  This  short 
duration  makes  the  frequency  content  of  the  pulse  on  the 

order  of  magnitude  of  1  MHz  which  is  in  the  region  where  the 

reactive  elements  of  the  cable  (capacitance  and  inductance) 
are  significant.  We  have  found  that  the  pulse  is  attenu¬ 
ated  from  lu  volts  to  about  5  volts.  The  triggering  thres¬ 
hold  of  the  ADS  is  3  volts,  thus  the  50%  attenuation  of  the 

signal  did  not  cause  any  problem.  An  additional  0.6  y  sec 
delay  was  observed,  corresponding  to  an  insignificant  rotor 
rotation,  and  we  concluded  that  the  triggering  process  was 


reliable. 


3.3  Data  Acquisition  Method 


Prior  to  starting  the  data  acquisition,  the  flight 
conditions  are  set.  This  was  done  by  adjusting  the  collec¬ 
tive  pitch  angle,  incidence  angle,  tip  speed  and  wind  velo¬ 
city.  Then  the  compensators  were  adjusted  according  to  the 
procedure  outlined  in  APPENDIX  A  to  yield  a  flat  frequency 
response  of  the  hot  film  and  data  output  system  up  to  about 
40u0  Hz.  The  next  step  was  to  statically  calibrate  the  hot 
film  anemometer  (outlined  in  sections  3./. 2. 2  to  3. 2. 2. 4. 2). 
After  the  calibration  constants  were  obtained,  they  were 
entered  manually  in  the  programs  (Appendix  B)  and  the  data 
acquisition  was  started  according  to  the  type  of  data  of 
interest,  i.e.  average  velocity  using  the  AVERAGE  program 
or  instantaneous  data  using  the  INSTVtL  program. 

The  main  purpose  of  this  work  was  to  develop  a  method 
of  measuring  instantaneous  velocities  in  the  wake  of  a 
helicopter  rotor.  However,  once  the  computer  and  the.ADS 
were  connected  to  the  wind  tunnel  we  used  it  for  average  data 
collection  as  well.  Both  programs  are  similar  and  the  only 
differences  are  in  the  averaging  and  storage  procedures. 

3.3.1  Average  Data  Collection 

In  this  mode  of  operation  the  data  was 
averaged  over  several  rotor  revolutions,  thus  the  angle 
azimuth  of  the  rotor  was  not  a  variable.  The  ADS  trigger¬ 
ing  was  done  directly  at  the  converter  by  an  exact  5  volt 
square  wave  obtained  from  a  TTL  output  of  a  WAVETECK  osci- 


lator. 


Since  the  estimated  significant  frequency  content 
was  4000  Hz  and  a  correct  reconstruction  of  a  wave  required 
at  least  two  samples  of  the  highest  frequency,  the  sampling 
frequency  was  selected  to  be  8000  Hz.  The  dimension  of  the 
array  BUFFER,  was  8192  (standard  dimensions,  see  Table  3.2), 
thus  the  array  was  filled  within  one  second,  which  corresponds 
to  about  40  revolutions  for  a  tip  speed  of  300  ft/sec. 

The  data  collection  was  done  in  a  predetermined  syste¬ 
matic  method,  which  is  basically  a  three  nested  do  loops, 
each  for  one  coordinate.  The  x  is  changing  the  fastest, 
then  z,  and  last  to  change  is  y.  Each  case  started  at 
x  =  28,  y  =  20,  z  =  12  and  from  this  point  on,  the  computer 
did  the  bookkeeping  and  instructed  the  operator  how  to 
proceed.  The  process  is  fully  continuous  from  previous 
days  (the  data  acquisition  for  one  case  might  take  about  a 
week)  with  no  information  needed  from  the  operator  except 
defining  the  status  of  the  new  day,  i.e.  is  the  present  data 
acquisition  a  continuation  of  previous  session  or  not?  If 
the  answer  is  yes,  the  computer  instructs  the  operator  where 
to  locate  the  probe,  and  if  the  answer  is  no,  the  computer 
instructs  the  operator  to  locate  the  probe  at  the  initial 
point.  Once  the  data  collection  starts  the  computer  will 
tell  the  operator  when  to  change  holders  and  warn  him  when 
the  probe  gets  into  the  area  in  which  the  rotor  is  too 
close  to  the  sensors. 
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At  each  point  in  the  wake,  the  data  was  averaged  (mean 
values)  for  each  channel,  the  effective  velocities  were 
calculated  and  transformed  into  wind  tunnel  coordinates. 

The  three  components  of  velocity  Vx ,  V  ,  Vz)  and  the  loca¬ 
tion  of  the  probe  (x-y-z)  were  stored  temporarily  on  the 
computer's  hard  discs.  At  the  end  of  the  day  the  data  col¬ 
lected  was  stored  on  a  floppy  disc  for  safety  reasons  so 
that  the  data  might  not  be  accidently  erased  and  lost.  After 
the  whole  case  was  completed,  the  file  was  transferred  to 
the  large  AMDAHL  470  computer  where  SPEAKEASY  language  was 
used  for  vector  plots  (Appendix  B) . 

3.3.2  Instantaneous  Data  Collection 

Using  the  assumption  that  the  wake  is  peri¬ 
odic  for  each  and  every  rotor  revolution  the  data  collected 
was  actually  the  average  instantaneous  velocity,  which  was 
averaged  at  every  five  degrees  increment  for  14  revolutions. 
For  this  case  the  ADS  had  to  be  triggered  each  time  the 
blade  had  passed  five  degrees  of  rotor  revolution.  At  a 
tip  speed  of  300  ft/sec  the  sampling  frequency  is  about  2750 
Hz  which  is  not  a  sufficient  sampling  frequency  for  a 
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complete  wake  reproduction.  Thus,  we  have  to  keep  in  mind 
that  the  values  measured  are  discrete  values  at  five 
degrees  increment  of  rotor  rotation.  The  five  degrees 
increment  triggering  signal  was  supplied  by  an  interface 
designed  and  built  for  this  project.  A  block  diagram  and 
elaboration  of  the  mode  of  operation  of  this  interface  are 
shown  in  Appendix  C. 

For  the  average  instantaneous  data  collected  two  methods 
were  developed.  1)  Computer  in  real  time  with  the  ADS. 

2)  The  analog  voltages  and  the  triggering  signal  recorded 
on  a  FM  analog  recorder  and  then  reproduced  at  the  ADS. 

The  second  method  was  developed  as  an  alternative  to  the 
first  one.  The  method  is  operational  but  has  not  been  used 
for  actual  data  acquisition. 

3. 3. 2.1  Computer  Mode  Operation 

The  inputs  to  the  interface  in  this 
mode  of  operation  were: 

1.  1/REV  pulse  from  a  magnetic  proximity  pickup  indicating 
\p  =  0. 

2.  72/REV  pulses  from  a  second  magnetic  proximity  pickup 
indicating  five  degrees  increments  in  rotor  rotation. 

The  output  from  the  interface  was  one  pulse  at  each  five 
degree  increments  for  15  revolutions.  (Figure  3.13) 

The  ADS  was  initiated  by  the  software,  but  the  digi¬ 
tizing  process  started  only  after  the  operator  signaled, 


TRIGGERING 

SIGNAL 


Figure  3-13  Triggering  Sequence  in  Computer  Mode 
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from  the  interface,  that  the  probe  was  at  the  correct  loca¬ 
tion.  Then,  the  ADS  was  triggered  every  five  degrees  as 
shown  in  Figure  3.13. 

After  the  digitizing  process  was  completed,  the  compu¬ 
ter  averaged  the  data  for  every  five  degree  increments  for 
14  revolutions.  The  three  average  voltages  were  transformed 
into  wind  tunnel  coordinates  in  the  usual  manner.  These 
velocities  and  the  probe  location  were  temporarily  stored 
in  an  encoded  way,  (see  Appendix  B)  on  the  computer's  hard 
discs,  until  the  case  was  completed.  The  data  was  then 
stored  permanently  on  floppy  discs  for  further  use  on  a 
reference . 


3. 3. 2. 2  Tape  Mode  Operation 

At  this  mode  of  operation,  the  inter¬ 
face  had  a  dual  purpose.  When  recording  the  data,  it  pro¬ 
duced  the  triggering  signals  for  each  five  degree  increments 
and  also  operated  the  FM  tape  recorder.  The  procedure  is 
as  follows:  the  operator  signaled,  from  the  interface, 
that  the  probe  was  at  the  correct  location,  then  the  device 
started  the  tape  recorder  and  after  0.5  seconds  delay  (to 
let  the  tape  get  to  its  steady  state  velocity)  the  analog 
data,  stepped  down  to  tape  recorder  maximum  allowable  input, 
was  recorded  on  three  channels.  On  the  fourth  channel,  a 
triggering  signal  was  recorded.  After  15  cycles  had  been  counted 
the  interface  stopped  the  recording  and  the  tape  was  shut 
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Figure  3-1  Recording  the  Triggering  Signals  in  Tape  Mode 
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Figure  3-14  Triggering  Sequence  in  Tape  Mode 
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down  (Figure  3.14a).  In  the  reproduction  process  the  three 
analog  channels  were  connected  to  the  ADS  and  the  fourth 
channel,  with  the  triggering  signals,  was  connected  to  the 
interface  whose  role  was  the  same  as  in  the  computer  mode 
(Figure  3. 14b) . 

After  the  digitizing  process  was  completed  the  data 
reduction  method  is  similar  to  the  procedure  outlined  in 
the  computer  mode. 


3.4  Data  Presentation 

The  numerical  values  of  the  velocities  and  their 
location  were  stored  on  floppy  discs  and  could  be  displayed 
on  any  terminal  or  printed  out  by  the  PR1NTR0NIX  high 
speed  printer.  However,  the  vast  amount  of  "numbers" 
associated  with  each  case  (1188  points  x  (72  x  3  numbers 
per  point)  =  256608  "numbers"  only  for  velocity  presenta¬ 
tion)  mode  the  use  of  data  visualization  a  necessity. 

For  data  visualization  the  TEKTRONIX  4014-1  graphics 
terminal  was  used.  This  terminal  has  a  refresh  mode  of 
operation  and  animated  pictures  of  any  plane  (xy,  yz,  zx) 
as  the  rotor  passes  by  were  created.  In  addition,  graph 
of  velocity  vs.  rotor  position,  or  velocity  vs.  probe 
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location  could  be  obtained  by  using  the  PLOT-10  plotting 
routines  available  in  the  PDP's  library. 

3.5  Measurement  Error  Analysis 

This  error  was  evaluated  experimentally,  in  two 
stages.  First,  the  velocity  calibrator  was  used  to  evalu¬ 
ate  the  velocity  and  angular  errors  because  of  imperfection 
of  the  sensors,  disturbances  due  to  the  prongs,  errors  in 
the  A/D  conversion  and  errors  of  the  fitted  calibration 
curves.  The  second  error  evaluation  was  done  in  the  wind 
tunnel,  to  evaluate  the  additional  errors  because  of 
misalignment  of  the  probe  in  the  wind  tunnel  and  because 
of  the  vibration  of  the  probe. 

3.5.1  Error  Evaluation  on  the  Velocity  Calibrator 
The  probe  was  installed  on  the  probe  holder 
in  the  calibration  nozzle  rig  so  that  all  the  sensors  were 
inside  the  core  of  uniform  velocity.  Then  the  velocity  of 
the  jet  was  set  to  velocity  values  between  10  and  70  ft/sec, 
and  the  probe  was  rotated  as  to  achieve  different  combina¬ 
tions  of  ^ .  The  velocity  was  varied  in  steps  of  10  ft/sec 
and  combination  of  the  directional  angles  were  set  at  each 
velocity.  The  anemometer's  analog  voltages  were  connected 
to  the  ADS  in  the  usual  way,  and  the  velocities  were  read 
from  the  terminal.  These  velocities  were  then  compared 
with  the  actual  velocities  determined  using  the  velocity 
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calibrator's  calibration  curve .  The  results  showed  a  maxi¬ 
mum  velocity  error  of  1  ft/sec  and  an  angular  error  of  +_  1°. 

3.5.2  Error  Evaluation  in  the  Wind  Tunnel 

For  this  case  the  probe  was  installed  in 
the  wind  tunnel  on  its  holders.  The  wind  velocity  was 
varied  in  the  same  range  as  on  the  velocity  calibrator. 
Vertical  and  horizontal  traverses  were  run  and  the  veloci¬ 
ties  obtained  were  compared  with  the  stream  velocity.  The 
magnitude  of  the  measured  velocity  was  within  1%  of  the 
stream  velocity.  However,  the  measured  stream  vector  had 
a  constant  5°  error  relative  to  the  X  axis  for  all  the 
points  measured  and  for  the  whole  velocity  range.  This 
error  was  due  to  imperfections  in  the  traverse  mechanism  and 
the  holders.  This  error  was  corrected  for  in  the  transfor¬ 
mation  equation  and  after  that  the  velocities  measured  in 
the  wind  tunnel  had  the  same  order  of  magnitude  of  error 
as  on  the  velocity  calibrator. 

An  extensive  investigation  done  by  the  author  for  an 
undergraduate  laboratory  design  course  (ME  581)  revealed 
that  the  velocity  error  due  to  vibration  at  the  sensor  is 


less  than  II. 


CHAPTER  IV 


TEST  RESULTS  AND  THEIR  PRESENTATION 

The  final  stage  of  this  work  was  a  preliminary 
measurement  of  a  z-y  plane  of  the  instantaneous  velocities 
in  the  wake  of  the  model  rotor  for  one  flight  conditions. 

The  purpose  of  this  stage  was  to  test  the  instrumentation 
and  mainly  to  develop  an  efficient  data  presentation  method, 
which  turned  out  to  tr  a  formidable  task.  The  parameters  of 
the  flight  condition  were  as  follows: 

Blade  Tip  Speed  ft/sec  300 

Advance  Ratio  0.06 

Blade  Twist,  degrees  8  (from  12.08  to  100%  R) 

Collective  Pitch,  (“5%  R)  8 

Rotor  Shaft  Tilt  Angle,  degrees  8 

The  instantaneous  velocities  were  measured  for  the  z-y 
plane  at  X/R  =  -  1.07.  This  plane  contains  66  points  and 
we  measured  the  three  components  of  velocity  every  five 
degrees  of  rotor  rotation  for  each  point.  Thus,  the  total 
number  of  velocity  information  for  the  entire  plane  was 
142S6.  (For  comparison  a  whole  case  of  average  velocity 
measurement  contained  "only"  3564  velocity  points.) 
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Therefore  it  is  clear  that  the  task  of  data  presentation 
was  not  a  simple  one. 

Our  first  inclination  was  to  use  the  vector  plots 
developed  for  the  average  data  presentation  (APPENDIX  B) 
for  each  azimuth  angle.  However  these  vector  plots  did  not 
transfer  to  the  viewer  any  information  typical  to  instan¬ 
taneous  velocities  that  is,  high  peaks  of  velocity  where 
vortices  were  encountered  by  the  sensors.  Figure  4-1 
shows  the  vector  plots  every  20  degrees  between  0  and  180 
degrees  azimuth.  The  plots  are  very  similar  and  very 
little  vector  movements  are  apparent  between  each  azimuth 
angle . 

Since  the  vector  plots  failed  to  transfer  the  desired 
information,  we  decided  that  a  presentation  of  the  indivi¬ 
dual  components  of  velocity  at  some  points  of  interest  was 
necessary.  Thus  the  plots  in  Fig.  4-2  and  4-3  were  created. 
These  plots  show  the  normal  induced  velocity  (V  )  as  a 
function  of  the  perpendicular  coordinate  to  the  rotor  at  an 
azimuth  angle.  For  this  presentation  a  complete  y  traverse 
was  made  for  two  azimuth  angles,  0  and  90  degrees.  The 
plots  showed  a  distinct  variation  of  the  velocity  with  y 
location  and  also  significant  changes  in  the  flow  pattern 
were  observed  between  the  plots  of  ip  =  0  and  -  90. 

The  third  type  of  plots  we  used,  were  graphs  showing 
the  velocity  variation,  at  a  point  in  the  wake,  as  a  func¬ 
tion  of  the  azimuthal  position  of  the  rotor.  These  graphs 
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are  shown  in  Figures  4-4,  4-5,  4-6  and  are  transferring 
the  most  information  about  the  flow. 
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CHAPTER  V 


DISCUSSION  OF  RESULTS 

The  plots  presented  in  Chapter  IV  were  a  preliminary 
attempt  to  check  the  instrumentation  sensitivity  and  to 
develop  data  presentation  methods.  Therefore,  the  discus¬ 
sion  in  this  chapter  is  limited  to  the  plane  measured 
(X/R  =  -1.07)  and  to  the  graphs  displayed. 

The  vector  plots  shown  in  Fig.  4-1,  displayed  flow 
characteristics  similar  to  the  average  velocities  measured 
in  Ref.  10,  that  the  roll  up  of  the  vortex  sheet  is  similar 
to  the  roll  up  of  the  wake  of  a  low  aspect  ratio  fixed 
win^.  Very  little  vector  movements  are  apparent  between 
each  azimuth  position  presented.  This  may  probably  be 
attributed  to  the  fact  that  the  tip  vortices  are  rotating 
in  the  X-Z  plane  thus  making  the  velocity  component  in  the 
Y  direction  relatively  small  and  its  contribution,  to  the 
velocity  vector  in  the  Z-Y  plane,  negligible. 

Figures  4-2  and  4-3  are  a  Y  traverse  at  two  azimuth 
angles,  \p  =  0  and  ijj  =  90  respectively.  The  flow  pattern  is 
similar,  i.e.  at  both  azimuth  angles  the  upwash  and  down- 
wash  occur  at  the  same  locations.  However  it  was  observed 
that  the  magnitudes  of  velocity  are  larger  on  the  advancing 

side  at  ip  =  90,  and  on  the  retreating  side  at  ip  =  0. 
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Figures  4-4,  4-5  and  4-6  are  the  time  histories  of  the 


three  components  of  velocity  related  to  blade  azimuth 
position.  Table  4-1  can  be  used  as  reference  for  the 
figures . 


POINT  # 

X/R 

Y/R 

Z/R 

p025 

1 

o 

0.27 

0.8 

p026 

-1.07 

0.27 

0.6 

p02  7 

-1.07 

0.17 

0.4 

p028 

-1.07 

0.27 

0.2 

p029 

-1.07 

0.27 

0.0 

p0  30 

-1.07 

0.27 

-0.2 

The  graphs  show  a  distinct  variation  of  the  velocity 
profile  with  azimuth  and  normal  distance  from  the  rotor. 

It  is  clearly  observed  that  12”  (Z/R  =  0.8)  below  the 
rotor's  tip  there  is  no  vortex  activity  indicating  that 
this  point  CP 02  5)  is  in  the  region  of  uniform  flow  and 
outside  the  rotor's  wake  boundaries. 

Fig.  4-4  at  Z/R  =  0  and  Z/R  =  -.2  show  the  distinct 
2/REV  pulses  expected  close  to  the  tip  of  a  two  bladed 
rotor.  At  these  points  the  measured  velocities  are  influ¬ 
enced  mainly  by  the  closeness  to  the  tip  vortices.  As  the 
measured  points  are  further  below  the  rotor  there  is  now  a 
combined  influence  of  the  bound  vortices,  the  shed  tip 
vortices  and  the  vortex  sheet,  creating  velocity  time 
histories  quite  different  from  the  2/REV  measured  closer  to 
the  blade  tip. 
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The  magnitudes  of  the  velocity  components  V  ,  V  ,  and 

x  y 

VT  as  observed  at  Z/R  =  0.0  and  Z/R  =  -0.2  indicate  that  the 
tip  vortex  rotates  in  the  X-Z  plane.  (V  and  are  much 
larger  than  V  ).  The  change  in  magnitude  of  V  in  Fig.  4-4 

y  * 

when  traversing  the  probe  from  Z/R  =  0.0  to  Z/R  =  -.2  indi¬ 
cates  that  the  vortex  rotates  in  a  counter  clockwise  direc¬ 
tion  and  that  the  Z  coordinate  of  the  vortex  core  lies 
between  Z  =  0  (p029)  and  Z  =  -3  CP 030)  .  The  reason  is  that 
at  Z  =  0  the  velocity  induced  by  the  tip  vortex  adds  to  the 
stream  velocity  and  at  Z  =  -3  (p030)  the  velocities  subtract 
creating  a  CCW  rotational  direction.  This  conclusion  is 
substantiated  by  the  fact  that  the  tip  vortex  has  to  rotate 
in  a  CCW  direction,  from  the  lower  side  of  the  blade  (high 
pressure)  to  the  upper  side  of  the  blade  (low  pressure) . 

Using  the  CCW  direction  of  rotation  of  the  vortex,  as 
explained  above,  the  velocities'  magnitude  at  Z  =  0  (p029) 
and  Z  =  -3  (p030)  in  Fig.  4-6,  indicate  that  the  measure¬ 
ments  did  not  traverse  thru  the  core  of  the  vortex,  but 
rather  on  one  side  only.  The  X  coordinate  of  the  vortex 
core  lies  to  the  left  of  the  X  coordinate  of  the  measured 
plane  (X/R  =  -1.07) . 

An  additional  interesting  feature  that  the  data  revealed 
was  a  difference  in  tip  vortex  strength.  Assuming  that  the 
collective  pitch  and  twist  of  the  blade  were  the  same  for 
each  blade  (our  measurements  were  within  +_  0 . 5  deg.)  the 
differences  in  vortex  strength  could  be  attributed  to 
differences  in  plane  of  rotation  between  the  two  blades. 
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If  one  of  the  blades  is  in  a  lower  plane  than  the  other, 
the  vortices  shed  by  the  blades  will  not  be  at  the  same 
location  and  the  velocities  induced  by  each  vortex  at  a 
given  point  in  space  will  be  different,  such  that  the 
furthest  one  will  have  lower  velocities  induced. 


CHAPTER  VI 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  instrumentation  and  the  methods  developed  were 
found  to  be  a  useful  tool  for  measuring  the  instantaneous 
velocities  in  the  wake  of  the  model  rotor.  The  results  of 
the  preliminary  tests  show  that  the  instrumentation  and  the 
computer  were  able  to  handle  the  vast  amount  of  measure¬ 
ments,  and  to  detect  important  characteristics  of  the  flow 
field. 

It  is  recommended  to  use  for  data  acquisition  a 
graphics  terminal  (TEKTRONIX  4014)  and  to  make  the  plotting 
routines  on-line.  Such  on-line  graphics  capabilities  will 
enable  the  operator  to  concentrate  more  in  the  regions  of 
interest  and  to  take  less  data  in  areas  of  low  vortex 
activity . 

For  better  confidence  in  the  collected  data  it  is 
recommended  to  have  the  computer  control  the  RPM  of  the 

rotor  and  the  dynamic  pressure  in  the  wind  tunnel.  The 
data  acquisition  would  be  done  only  when  the  computer  sig¬ 
nals  that  the  conditions  in  the  wind  tunnel  are  correct. 
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Appendix  A 

A  Method  of  txtending  the  Frequency  Response  of  Hot 
Film  Anemometers  Beyond  Their  Natural  Capacity 

A.l  Determination  of  the  Significant  Frequency  Content  of 
the  Wake  of  a  Two  Bladed  Rotor. 

The  highest  frequency  that  the  measuring  system 
encounters,  will  be  when  the  vortex  sheet  or  a  tip  vortex 
passes  by  the  sensors.  For  this  case  the  velocity  has  an 
abrupt  change  in  magnitude,  similar  to  a  pulse,  and  a 
high  frequency  content  proportional  to  the  duration  of 
the  pulse.  The  shorter  the  pulse  is,  the  higher  will  be 
its  frequency  content  . 

In  order  to  get  an  estimation  of  what  significant 
frequencies  will  be  encountered  in  the  wake,  a  two  pulse 
per  revolution  created  by  the  tip  vortex  of  each  blade 
was  simulated  by  using  the  measurements  of  Ref.  32.  It 
was  assumed  that  the  tip  vortex  is  leaving  the  blade  at 
some  instant  and  is  "travelling”  with  the  wind  velocity. 
This  assumption  enabled  us  to  estimate  the  duration  of  the 
pulse  from  its  length  as  represented  in  Ref.  32.  The 
magnitude  of  the  pulse  does  not  have  any  effect  on  its 
frequency  content.  Fig.  A-l  shows  the  wave  shape  used  in 
the  Fourier  analysis  to  determine  the  frequency  content. 

It  was  assumed  also,  that  between  pulses  nothing  happens 
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Figure  A-l  Simulated  2/REV  Vortex  (Ref.  32) 


a.  100th  Harmonics  b.  50th  Harmonics  c,  25th  Harmonics 


Figure  A-2  Harmonics  Representation  of  the  Velocity  Pulse 
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in  the  wake.  This  assumption,  which  is  not  correct,  will 
yield  higher  harmonics. 

The  fundamental  frequency  is  the  frequency  of  the 
rotor,  which  is  10  Hz  in  the  case  of  Ref.  32,  and  each 
harmonic/  is  an  integer  multiplication  of  this  frequency. 

Using  a  canned  Fourier  program  it  was  found  that  using 
the  case  of  100th  harmonics  the  spikes  were  reproduced  with¬ 
in  99.9%,  the  50th  harmonics  creates  a  10%  error  only  at 
the  peak  and  the  location  of  the  peak  is  accurate.  The  25th 
harmonics  is  far  from  the  original  shape  but  it  is  still 
quite  a  good  approximation  to  the  velocity  pulse.  Fig.  A-2 
shows  the  positive  part  of  the  pulse  superimposed  on  the 
harmonics  representation. 

In  our  case  the  frequency  of  the  rotor  was  37.5  Hz.  It 
was  decided  to  strive  for  the  highest  harmonics  representa¬ 
tion  possible,  i.e.  100th  harmonics.  This  implied  that  the 
required  frequency  response  of  the  hot  film  anemometer  had 
to  be  flat  up  to  3750  Hz.  However,  the  0.002"  diameter  hot 
film  sensor  had  a  flat  frequency  response  up  to  about  1000 
Hz,  which  is  considered  high.  In  order  to  improve  the 
frequency  response  of  the  hot  film  anemometer  it  was  decid¬ 
ed  to  design  a  compensator  to  obtain  the  desired  range  of 
frequency.  Sections  A-2  and  A-3  outline  the  method  and 
procedures  developed. 
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A.  2  Direct  Determination  of  Hot-Film  Anemometer  Frequency 
Response 

The  method  presented  here  is  based  on  an  idea  from 
Reference  21.  Let  G(iw)  be  the  sinusoidal  transfer  function 
of  a  dynamic  system,  Q^(ia>)  the  input  spectrum  to  that 
system,  and  Q  (iu>)  its  output  spectrum.  Then  we  can  write 

Qo(ia0 

-  oft 1ST 

A  simpler  form  of  equation  (A-l)  is  one  using  decibel 
coordinates,  then  equation  A-l  becomes: 

201og10(G(iu))  =  201og(Qo(iu»))  -  201og CQ± (iw) )  (A-2) 

The  logic  behind  this  method  is  to  excite  the  system  with 
a  signal  of  known  frequency  spectra  to  measure  the  output 
spectra,  and  then  to  subtract  these  two  spectras.  An 
additional  practical  simplification  will  be,  if  the  input 
spectrum  is  constant,  up  to  a  maximum  frequency  (white 
noise),  then  the  output  spectrum  when  normalized  represents 
actually  the  sinusoidal  transfer  function  of  our  system. 

This  can  be  written  as 

201og1£)G(ia))  =  201og10Qo  (iw)  -  N  (A-3) 

e 

where  N  is  the  constant  input  spectrum.  In  decibals  N 
just  shifts  the  whole  output  spectrum  up  or  down. 

A  special  fan  was  designed  and  built  for  this  purpose. 
The  fan  is  a  4"  impeller  with  60  blades  attached  to  an 
electric^  motor  with  variable  speedf  up  to  10000  RPM. 
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The  impeller  produced  velocity  white  noise  within  +_  0,5  db 
up  to  lOKHz  at  10000  RPM.  The  spectrum  was  measured  by  a 
TSI  0.0004"  diameter  hot  wire  with  flat  frequency  response 
up  to  14  KHz  measured  by  an  indirect  method  discussed  in 
Ref.  30.  Figure  A-3  shows  the  impeller  and  the  reference 
hot  wire,  and  Figure  A-4  shows  the  input  spectra,  calculated 
by  a  Nicolet  660-A  digital  spectrum  analyzer. 

Once  the  white  noise  characteristics  of  the  impeller 
spectrum  was  established  we  proceeded  to  obtaining  the  hot 
film's  anemometer  break  point. 

The  tested  hot  film  sensor  was  located  at  the  same 
location  where  the  references  spectrum  was  obtained.  The 
measured  spectrum  was  the  hot  film's  frequency  response  as 
explained  before.  The  point  where  the  assymptote  of  the 
flat  frequency  response  meets  the  assymptote  to  the  roll 
off  indicates  the  hot  film's  break  point  (f^p) .  (Fig.  A-6a) 
A. 3  The  Design  of  a  Lead  Compensator  to  Extend  the  Frequen¬ 
cy  Response  of  the  Hot -Film  Anemometer 
Once  the  break  point  of  the  hot  film  anemometer  was 
established,  as  described  in  Section  A-2,  a  lead  compensa¬ 
tor  was  designed  to  have  a  roll  up  to  compensate  for  the 
hot  film's  roll  off. 

The  authors  of  Ref.  27  have  shown  that  the  hot  film 
rolls  off  with  a  slope  of  -30  db/DECADE  ignoring  reactive 
elements.  This  kind  of  slope  is  very  hard  to  obtain  by 
means  of  simple  operational  amplifier  circuits.  However, 


77 


the  authors  of  Ref.  27  have  mentioned  that  the  actual  slope 
measured  was  steeper.  Experiments  run  for  this  work,  using 
the  method  described  in  Section  A-2,  revealed  that  the  hot 
film  rolls  off  with  a  slope  of  -42  to  -47  db/DECADE.  This 
slope  is  closer  to  the  +40  db/DECADE  slope  obtainable  with 
OP-AMP  circuits.  A  typical  compensator  designed  and  built 
for  this  work  is  shown  in  Figure  A-5. 


Rx  =  10  t  250  KO  R2  =  15  KO 

C1  =  5110  pF  C2  *  90.9  pF 

OP  =  CD471  OPERATIONAL  AMPLIFIER 
Figure  A-5  Lead  Compensator  Circuit 


The  transfer  function  of  the  compensator  is: 


(A-3) 


The  numerator  is  a  second  order  system  with  damping  ratio 
of  one  and  a  roll  up  of  40  db/DECADE.  The  denominator  is 
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a  second  order  filter  needed  to  attenuate  the  high  frequency 
noise  associated  with  differentiators.  In  addition,  three 
antialising  Chebychev  filters  were  used  in  the  ADS  to  fur¬ 
ther  attenuate  the  signal  for  frequencies  above  5  KHz.  The 
break  point  of  the  compensator  is  adjustable  between  400  to 
4000  Hz,  according  to  the  measured  break  point  of  the  hot 
film  anemometer  system.  The  compensator's  adjustment  is  as 
follows:  the  hot-film's  break  point  (f^p)  is  obtained  as 

described  in  Section  A-2.  We  have  found  that  the  best 


compensation  is  obtained  when  the  damping  ratio  of  the  hot 
film  anemometer  system  is  adjusted  to  a  higher  value  than 
the  optimum  required  (s  =  0.7)  for  a  good  response.  This 
is  done  by  adjusting  the  offset  voltage  in  the  anemometer 
circuit.  The  higher  damping  ratio  causes  a  more  sluggish 
response  of  the  hot  film  but  if  the  compensator's  break 
point  (fc)  is  adjusted  to  approximately  1.5  f^p,  the  total 
response  is  flat  within  +_  0.5  db  up  to  5000  Hz.  From  the 


value  of  f  ,  Rp  can  now  be  calculated  by  means  of  equation 


R  =  -  -1-  .  1 

K1  27Tfc  Cx  2 IT  c  1  •  5  fHF) 


(A-4) 


where  Cp  »  5110  x  10*  F.  Any  adjustments  in  Rp  effects 
the  gain  of  the  compensator  as  shown  in  equation  A-3.  The 


value  of  R2  was  chosen  so  that  the  gain  changes  between  1:3 
to  3:1  at  the  extreme  designated  break  points  of  400  and 
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problems  where  as  too  low  a  value  saturates  the  OP-AMP, 

To  overcome  these  two  problems  the  potentiometer  Rq  was 
introduced  to  increase  or  decrease  the  signal  according  to 
the  existing  situation.  Figure  A-6  shows  the  response  of 
one  of  the  channels  to  a  velocity  white  noise,  a,  without 
compensation  and  b.  with  the  lead  compensator  present. 


<  23 


Figure  A-6b  Hot -Film  Anemometer  Response  With  Compensator 
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Appendix  B 

Computer  Programs  for  Instantaneous  and 
Average  Data  Acquisition. 

B.l  Instantaneous  Velocities  Acquisition  Program  (INSTVEL) 
INSTVEL  is  an  interactive  program  used  for  instantan¬ 
eous  velocities  data  acquisition.  ihe  program  controls 
the  analog  to  digital  converter  system  (.ADS)  by  means  of 
subroutines  SETAD1  and  GETAD.  The  ADS  sends  to  the 
computer  the  discrete  voltage  values,  as  measured  by  the 
3-D  hot  film  anemometer,  every  five  degrees  of  rotor 
revolution  and  for  15  revolutions.  The  program  selects 
the  five  degree  increment  of  each  revolution  from  the 
array  BUFFER,  averages  these  values  (using  the  subroutine 
AVERAG)  and  then  transforms  the  measured  voltages  into 
the  effective  cooling  velocities  using  the  calibration 
constants  of  each  channel.  In  the  next  step  the  program 
calculates  the  total  velocity  vector  (VEL)  and  its  compo¬ 
nents  in  the  probe  coordinate  system  iterating  between 
2 

VEL,  k  and  <{k  using  subroutine  UKITER  and  the  functions 
UKl  for  sensor  A,  UK2  for  sensor  B  and  UK3  for  sensor  C. 
Once  the  velocity  components  in  probe  coordinate  system 
were  obtained  they  are  transformed  into  wind  tunnel  coor¬ 
dinates.  The  velocities  so  obtained  are  encoded  into  a 
one  nine  digit  integer  number  and  stored  in  a  direct  access 
mode.  This  is  done  as  follows:  the  three  velocity  compo¬ 
nents  are  multiplied  by  10,  to  make  the  maximum  velocity  a 
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three  digit  value,  and  a  constant  value  of  499  is  added  to 
each  component  to  make  the  velocity  positive.  Now  IVX  is 
multiplied  by  10^,  IVY  by  103,  and  the  three  values  are 
added.  In  the  program  this  is  done  as  described  below: 

IVX  =  IFIX(Vx  *  10.  +  0.5)  +  499 

IVY  =  I FIX (V  *  10.  +  0.5)  +  499 

IVZ  =  IF  IX (V  *  10.  +  0.5)  +  499 

NVEL  =  IVX*106  +  IVY  *  103  +  IVZ 

Thus,  NVEL  is  a  nine  digit  integer  number  which  requires 
less  memory  space  than  the  three  real  values  of  V^,  V  ,  V^. 
In  addition,  the  direct  access  storage  mode  compressed  the 
data  even  more.  This  special  way  of  data  handling  enabled 
us  to  reduce  the  number  of  floppy  discs  required  to  store 
the  data  for  one  case,  from  6  to  one  floppy  disc.  The 
location  of  the  probe  (x-y-z)  was  coded  and  stored  in  a 
similar  way.  The  velocity  coding  and  decoding  caused  a 
maximum  error  of  0.05  ft/sec  which  was  well  within  our 
measurements  error. 

B.1.1  Averaging  Method  for  INSTVEL 

The  averaging  was  done  on  the  digitized  voltages 
in  order  to  reduce  the  number  of  data  points  required  for 
velocity  calculations.  The  computer  selected  all  the  data 
points  at  each  azimuthal  angle  (every  72nd  element  in  the 
array  BUFFER)  and  calculated  the  mean  of  these  point  using 
equation  B-l. 
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B.2  Average  Velocities  Acquisition  Program  (AVERAGE) 

The  AVERAGE  program  is  similar  in  many  ways  to  the 
INSTVEL  program.  Thus,  only  the  significant  differences 
are  elaborated  here. 

For  the  average  velocities  the  azimuth  angle  of  the 
blades  is  not  relevant.  Therefore,  the  averaging  was  done 
over  several  rotor  revolutions.  The  ADS  is  triggered 
continuously  by  a  TTL  output  from  a  Waveteck  Oscillator 
until  the  array  BUFFER  is  filled.  A  new  array  (DATA)  is 
defined  which  includes  all  the  digitized  value  of  one 
channel  at  the  time.  Then,  the  mean  is  obtained  for  each 
channel  and  the  values  are  processed  in  the  usual  manner. 

The  three  velocity  components  (Vx ,  V  ,  Vz)  and  their 
location  are  stored  in  a  FORMAT  (317,  3F7.2)  on  a  hard  disc. 
When  the  case  is  completed  the  whole  data  set  is  transferred 
to  a  magnetic  library  available  with  the  large  ADMAHL  470 
of  The  Ohio  State  University,  so  that  the  vector  plotting 
routines  (Section  B.3)  could  be  used. 

B.2.1  Averaging  Method  for  AVERAGE 

The  digitized  voltages  were  averaged  so  that 
memory  space  and  computational  time  could  be  saved.  For 
this  case  the  data  points  of  each  channel  were  averaged 
over  40  rotor  revolution,  which  corresponds  to  2048  data 
points.  The  program  selected  the  first  2048  elements  from 
the  array  BUFFER  for  each  of  the  channels  and  using  sub¬ 
routine  AVERAGE  the  mean  was  found  by  means  of  equation  B-2. 
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B.3  SPEAKEASY  Programs  for  Average  Velocities  Vector  Plots 
The  SPEAKEASY  language  proved  to  be  a  very  useful  and 
convenient  means  of  data  presentation.  The  following 
sections  shows  the  programs  and  the  vector  plots  obtained. 
The  only  difference  between  each  program  is  the  number  of 
vectors  used  for  the  plot  according  to  the  planes  plotted, 
and  the  indicator  (I)  used  to  select  the  proper  vectors. 
B.3.1  Vector  Plots  of  Z-X  planes. 

For  this  plane  I  is  defined  as: 

I  =  GRID(1 ,1188,18) 

and  66  vectors  were  used  for  each  plot.  Figure  B-l  shows 
these  planes  plotted,  and  the  additional  information 
printed,  by  the  computer  on  each  graph. 
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62  TEXTPUT<TXT2,10,3> 

63  TEXTPUTC20  FT/SEC"* 

64  1*1+1; PRINT  NAX< I > 
*66  FREE  TECK 

■X 


10,-6) 


[iE 


forumo  x/r 

orr 

forumro  x^» 

or 

__  20  ft /see 

v/n— i.r> 

_  e»  ft ^ see 

v/»*f.S33 

FORyaRO  x/» 


__  a*  r t^scc 


i.s  e.a  t.a  *0.8 

FOOUMD  X/R 

-I.S  “*.« 

•ft 

FORWARD  X/R 

M 

.  M  Ff/SfC 

V/ff— «.S?? 

__  e»  ft/scc 

V/R«l .#7 

rORuMD  X/R 

err 

FORUMO  X/R 

R 

_  eo  FT  'SCC 

V/R»-R.«*7 

_  »  ft/ sec 

v/R-l.M 

Figure  B-l 


Velocity  vectors  in  X-Z  plane 
for  y-0.06,  8_-%  ■  8°,  and 
i-8°.  Note  tnat  coordinates 
depict  probe  position  only. 
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B.3.2  Vector  Plots  of  X-Y  Planes. 

For  this  case  the  initial  I  is  defined  as 
I  =  GRID(1, 108,1) 

and  108  vectors  were  used  for  each  plot.  Figure  B-2 
these  plots. 
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EOITING  ARROMXY 

1  PROGRAM 

2  HODEPRINTCOFF) 

3  GET  CORTM 

4  XX=CORTM<  I  / ) 

3  FREE  CORTM 

6  AA*2/3 

7  Q«XX<1,3>/13 

8  Q*R0UN0ED<Q,3> 

9  Q=BASE<  Q# 10 ); Q=CHARACTER(  Q > 

10  TEK(120) 

11  X*-XX<  1 2 AA 

12  X1*X-<XX<,4H3.06)/6*AA 

13  MHERE( XX( #  4 ) . EQ . 0 )  X1«X 

14  Z=XX(,1>*AA 

15  Z1*Z+XX<  » 5  )/6$AA 

16  MHERE( XX< / 4 ) . LT . -90  >  Z1«Z 

17  B*.6;C*.7 

18  FREE  XX 

19  DZZlsABS< ABS(Z)'ABS(Z1 >) 

20  DXX1=ABS<ABS<X>-ABS<X1» 

21  X2«X*B*DXX1 

22  X3«X*C*0XX1 

23  Z2»Z*C*DZZ1 

24  Z3«Z*B*DZZ1 

23  MHERE(X.GT.Xl)  X2-X-6SDXX1 

26  HHERECX.CT.Xl)  X3*X-C*DXX1 

27  MHERECZ.GT.Zl >  Z2«Z-C*DZZ1 

28  MHERECZ.GT.Zl >  Z3*Z-BtDZZi 

29  HSCALE-C-33.45) 

30  USCALE-C -30,30) 
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31  JOINC XI, Z1 ,X2,Z2) 

32  J0INCX1.Z1.X3.Z3> 

33  JOINCX.Z.Xl.Zl) 

34  FREE  X  Z  XI  Z1  X2  Z2  X3  Z3 

35  FREE  0XX1  0ZZ1 

36  JO INC -10,-20,-10.30 ) 

37  JOINC-l0.-2O,8,-20> 

38  JOINC -10, 30,8,30) 

39  JOINC 8,30,8,-20 > 

40  XTIC*ARRAYC6: )-10 

41  XTIC1«XTIC*.5 

42  N=GRIDC0,5, 1 ) 

43  ZTIC1=ARRAYC6:>-16+12*N*AA 

44  JOINC XTIC,ZTIC1,XTIC1,ZTIC1 > 

45  N-GRIOC 0,4,1) 

46  ZT I C 1 « ARRA YC 5  * )-8+6*N*AA 

47  XTIC=ARRAYC5=  >-20;XTICt*XTIO.5 

48  JOINC ZTIC1,XTIC,ZTIC1.XTIC1 ) 

49  JOINC 10,-16, 10, -16+20/6*AA> 

58  D*-12iE*-24 

51  TEXTPUTC "1 .6*, 0,-17) 

52  TEXTPUTC "0.8", 0,-9) 

53  TEXTPUTC "0. 0".D. -I > 

54  TEXTPUTC "-0.8", 0.6) 

55  TEXTPUTC "-1. 6", 0,14) 

56  TEXTPUTC "-2. 4", 0,22) 

57  TEXTPUTC "FORHARO" , 0-2,  -18 . 67 > 

58  TEXTPUTC "X/R". 0-2,0) 

59  TEXTPUTC "AFT ",0-2, 24) 

60  TXT10=ARRAYC3, 1 « "Z/R" >i TEXTPUTC TXT10,0,E-2> 

61  TEXTPUTC "20  FT/SEC", 10,-11 > 
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62  TEXTPUTC  "-8 . 8" , -9, E-l >j TEXTPUT<  "8 . 8* , -1 , E > 

S  tIxtpI!t!  Wxlf-' >,TEXTPUT<  ■••4,'5'E> 

65  TXT12*"Y/R*",Q 

66  TEXTPUTC  TXT12#  18/ 16  > 

6?  1*1+188 

68  PRINT  HAX<I) 

$69  FREE  TECK 


29  FT/SCC 


2/9*9, 2 


_ 29  FT/SEC  2/9*9 


A 


_ 29  FT*- SEC  .  Z/R— 9.2 


A 


Figu 

Velocity  ve 
for  V«0.06, 
i«8° .  Hot 
depict  prob 


a 


a 


_  29  FT/SEC 


2/R— 9.4 


A 


_ 29  F Tv SCC  Z/R—9.6 


A 


_  29  FT/SEC  2/R— 9.2 


A 


re  B-2 

ctors  in  X-Y  plane 

®75%8'  "  ®°'  *n<5 
e  tnat  coordinates 

e  position  only. 
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B.3.3  Vector  Plots  of  Z-Y  Planes. 

For  these  planes  the  indicator  I  was  more  complex 
since  the  needed  points  were  the  first  18  points  then  108 
had  to  be  shipped  and  so  on,  therefore  the  simple  GRID 
statement  could  not  be  used.  A  small  computer  program 
was  written  to  create  the  I  for  this  case,  and  it  was 
stored  on  our  magnetic  library  for  future  use. 

Figure  B-3  shows  the  Z-Y  planes  obtained  by  the  plot¬ 
ting  program. 


EDITING  ARROHZY 
1.0  PROGRAM 
2.0  HQDEPRINT< OFF) 

3.0  AA*2/3 
4.0  GET  CORTH 
5.0  XX*CORTM<  I » ) 

6.0  FREE  CORTH 
7.0  Q*XX(1,2> 

8.0  Q*-Q/15 
9.0  Q*ROUNDED( Q*  3  ) 

10.0  Q=BASE(Q* 10) 

11.0  QSCHARACTER( Q ) 

12.0  X-XX< / 3 )*AA 

13.0  Xl«X+< XX<  /6)-5.57 V6TAA 

13.5  UHEREC  XX( *  6  ) . EQ . 0 )  X1*X 

15.0  Z*XX< , 1 )SAA 

16.0  Zl«Z+XX(,5)/6tAA 

18.0  B=.6;C=.7 

19.0  FREE  XX 

20.0  DZ21*ABS( ABS<  Z )-ABS< Z1 ) ) 
21.0  DXX1*ABS<ABS<X>-ABS<X1>> 
22.0  X2=X+B*DXX1 
23.0  X3=X+C*DXX1 
24.0  Z2*Z+C*0ZZ1 
25.0  Z3sZ+BtDZZl 
26.0  MHERE(X.GT.X1 )  X2«X-B*DXX1 
27.0  HHERE(X.GT.Xl)  X3*X-C*DXX1 
28.0  HHERECZ.GT.Zl >  Z2-Z-C*0ZZ1 
29.0  HHERE(Z.GT.Zl)  Z3«Z-BSDZZ1 
30.0  HSCALE-< -35,43) 

■*. 


31  USCALE»< -30  , 30  > 

32  TEKC120) 

33  J0IN(X1,Z1,X2,Z2> 

34  J0IN<X1,Z1,X3,Z3> 

35  J0IN(X,Z,XI,Z1> 

36  FREE  X,Z,X1,21,X2,X3,Z2,Z3 

37  J0IN<-16. 67, -20,-16. 67, 30) 

38  JOINC-16. 67, -20, 16.67,-20) 

39  JOIN( -16.67,30, 16.67,30) 

48  JOIN< 16.67,30, 16.67,-20) 

41  XTIC=ARRAY<6: )-16.67 

42  XTIC1=XTI0.5 

43  N=GRID(0,5, 1 ) 

44  ZTIC1*ARRAY(6: )-16+12*N*AA 

45  J0IN(XTIC,ZTIC1,XT1C1,2T1C1 ) 

46  XTIC*ARRAY(6!  )-20;XTlCl*XTIC%5 

47  ZTIC1=ARRAY<6- )-16M2*N*AA 

48  JOIH< ZTIC1,XTIC,ZTIC1,XTIC1 > 

49  JO I N<  20 , - 1 6 , 20 , - 1 6-20^6 tAA ) 

50  D=-18;E*-24 

51  TEXTPUK "1 .6%D,-17) 

52  TEXTPUK  "0. 8% D, -9) 

53  TEXTPUK  "0.0",D,-1 ) 

54  TEXTPUK --0.8*, 0,6) 

55  TEXTPUK  "-1. 6% 0,14) 

56  TEXTPUK *-2. 4", D, 22) 

57  TXT7* "FORWARD ■  ;  TEXTPUK  TXT7 ,0-2,-18.67) 

58  TXT8*-X/R- ; TEXTPUK  TXT8,0-2, 0  > 

59  TXT9* " AFT “ ; TEXTPUT<  TXT9 , 0-2 , 24  > 

60  TXT 1 0* ARR A Y<  3 , 1  *  ■  Y /R 1 ) ;  TEXTPUK  TXT  1 0 , 0 ,  E-2  > 

61  TXT11--20  FT/SEC- ; TEXTPUK TXT1 1,20, -15) 
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62  TEXTPUT< - 1 . 6" * -16, E ) ) TEXTPUT< *8.8*,-8,E) 

63  TEXTPUK  "8 . 8' ,  8 ,  E ) ;  TEXTPUK " -8 . 8" , 8 , E- 1 > 

64  TEXTPUK  M.6MS,E-1) 

65  TXT 1 2*ARRAY< 5 , 1  * “ AOUAN " > 

66  TEXTPUK  TXT12, 18,  E-2) 

67  TXT 1 3=ARRAY<  4,1s "RETR ■ > 

68  TEXTPUK TXT13, -13, E-2) 

69  TXT14=-Z/R=',0;TEXTPUT<TXT14,28,16) 

70  I=i+18;PRINT  HAX< I ) 

*71  FREE  TECK 

*4 


Appendix  C 

ADS  and  Tape  Recorder  Triggering  Circuit  -  Interface 

The  purpose  in  designing  and  building  this  device, 
was  to  provide  a  simple  method  of  measuring  the  rotor's 
azimuthal  angle  at  five  degree  increments  and  to  trigger 
the  ADS  at  each  increment  for  a  predetermined  number  of 
revolutions.  An  additional  mode  of  operation  wad  designed 
also,  that  can  be  used  when  the  analog  data  is  recorded 
first  on  a  FM  analog  tape  recorder. 

As  mentioned,  there  are  two  modes  of  data  acquisition: 
1.  Computer  Mode,  2.  Tape  Mode.  The  modes  of  operation 
are  selected  by  the  switch  SW  as  shown  in  Figure  C-l. 

Prior  to  starting  data  acquisition  the  number  of  revolutions 
(N)  is  set,  for  the  case  presented  in  this  work,  N  was  15. 

In  the  computer  mode  the  inputs  to  the  interface  are 
a  1/R  and  a  72/R  signals  from  two  proximity  magnetic  pick¬ 
ups.  The  comparators,  Cl  and  C3  are  shaping  the  approxi¬ 
mately  sinusoidal  pulses  from  the  magnetic  pickups,  into 
positive  square  pulses  whose  duration  is  the  duration  of 
the  pulses  from  the  pickups.  The  role  of  the  monostables, 
MSI  and  MS2 ,  is  to  create  a  pulse  of  constant  duration  and 
height  whenever  there  is  an  output  pulse  from  the  compara¬ 
tors.  The  operator  signals  that  the  digitizing  process  can 
start  by  pushing  the  push  button  PB.  This  signal  energizes 
the  flip  flop  (FF)  FF2  which  sends  a  constant  binary  1  to 
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the  nand  gates  (N)  N1  and  N2.  N1  energizes,  initially,  FF1 
and  also  sends  the  1/R  signal  to  the  counter  which  compares 
the  number  of  signals  coming  from  N1  with  the  preset  value 
of  N  revolutions.  FF1  remains  energized  as  long  as  the 
number  of  signals  from  N1  is  less  or  equal  to  N.  Since  FF1 
sends  a  constant  binary^  to  N3,  it  will  transfer  all  the 
72/R  pulses  to  MS8.  MS8  is  energizing  the  drive  device 
which  creates  a  10  volt  signal  of  1  usee  duration  each  time 
there  is  a  72/REV  signal.  After  15  revolutions  have  been 
counted  the  process  stops,  all  the  FF's  are  reset,  and  the 
system  awaits  the  new  initiating  signal  from  PB. 

In  the  tape  mode  the  interface  has  a  dual  purpose,  one 
in  the  record  mode  and  the  second  in  the  reproduction  mode. 

In  the  record  mode  SW  is  in  the  same  position  as  is 
the  computer  mode  and  the  inputs  1/R  and  72/R  are  the  same. 
The  A/D  triggering  is  disconnected  and  one  channel  of  a  FM 
analog  tape  recorder  is  connected  to  tape  out.  The  process 
is  the  interface  is  similar  to  the  computer  mode  operation, 
except  that  now,  the  signals  from  N1  and  N2  are  flowing  to 
MS8  and  MS9  respectively.  Then  MS8  creates,  thru  G1 ,  a  10V 
pulse  representing  a  1/REV  signal,  and  MS  (creates,  thru 
G2) ,  a  5V  pulse  representing  a  5  degree  increment  in  rotor 
rotation.  In  addition  the  operation  of  the  tape  is  con¬ 
trolled  by  the  interface.  After  the  operator  has  pushed 
PB  the  tape  is  operated  forward  and  after  0.5  sec  delay 


Figure  C-l  Interface  Generalized  Block  Diagram 
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the  recording  starts.  The  reset  signal  from  the  COUNTER 
energizes  MS7  which  stops  the  tape  completely. 

In  the  reproduction  process  the  fourth  channel  of  the 
tape  recorder  is  connected  to  TAPE  IN  and  SW  is  switched 
to  tape  mode.  C2  is  adjusted  to  react  only  to  the  10V 
signals  and  C4  reacts  to  the  5V  signals.  The  A/D  output 
of  the  interface  is  connected  to  the  triggering  of  the  ADS. 
When  the  operator  pushes  PB  the  same  process  occurs  as 
explained  in  the  computer  mode.  In  addition  the  tape  is 
operated  in  the  manner  explained  in  the  recording  mode. 

Figure  C-l  is  a  generalized  block  diagram  of  the  inter¬ 
face,  and  does  not  include  fine  details  of  the  design  which 
were  needed  to  make  the  system  work. 


